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■ Abstract. In this contribution, I will give a brief survey of present techniques to treat the bound 

state problem in relativistic quantum field theories. In particular, I will discuss the Bethe-Salpeter 
equation, various quasi-potential equations, the Feynman-Schwinger representation, and similarity 
transformation methods for Hamiltonian approaches in light-front quantization. Finally, I will 

\ comment on a related similarity transformation in the usual equal-time quantized theory. 



One of the most striking and ubiquitous phenomena in nature is the binding of matter. 
It is found on essentially all length scales, from galaxies to quarks, and is even speculated 
to be crucial for the description of physics much below the length scale of hadronic 
physics. We currently have a good understanding of the binding mechanism in crystals, 
molecules, and atoms. The non-relativistic Schrodinger equation provides a satisfactory 
description at these scales, and small corrections due to relativistic and virtual particle 
creation effects can accurately be treated in perturbation theory. However, at smaller 
length scales relativistic effects become increasingly important, and the description 
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through the non-relativistic Schrodinger equation is not accurate enough. It is hence 
inevitable to employ quantum field theory (QFT), the best currently known description 
of microscopic physics. 

In the following, I will give a brief survey of the most popular equations which have 
been employed for bound state calculations in relativistic QFT. Attention will be re- 
stricted to equations which are rooted in field theory, in particular, I will leave out the 
generalizations of relativistic one-particle quantum mechanics to two or more particles, 
as well as equations with phenomenological input. Equations that are derived from field 
theory can naturally be divided into two classes, those based on the manifestly covariant 
Lagrangian, path-intregal formulation of QFT, and those built on a Hamiltonian formu- 
lation in a perturbatively defined Fock space. 

In the Lagrangian approach to QFT, which is particularly convenient for the calcula- 
tion of transition amplitudes, a bound state of two constituents appears as a pole in the 
2-particle scattering amplitude, or of the field-theoretic 4-point correlation function. 
The pole is interpreted as the pole in the propagator of the bound state as illustrated in 
Fig. 1, and its position determines the bound state mass (via the usual pole mass defini- 
tion). From the figure it is clear that near the pole the 4-point function has the form 
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FIGURE 1. The contribution of an AB-bound state to the 4-point function for the scattering of particles 
A and B. 



where P = pa + Pb = p'a + Pb * s me total momentum, and finite contributions have been 
suppressed. The function 



is the so-called Bethe-Salpeter wave function and encodes the structure of the bound 
state | yr). 

It is a well-known fact that no contribution of finite order in perturbation theory, 
or, diagrammatically speaking, no finite number of Feynman diagrams, can ever give 
a bound state pole in the 4-point function. Hence it is imperative in the Lagrangian 
approach to sum up an infinite subset of diagrams. It is this property that makes the 
calculation of bound states a non-perturbative problem in QFT, and the choice of 
an appropriate subset of diagrams and its summation present formidable problems in 
practice. 

The first equation that was devised to cope with this problem is the Bethe-Salpeter 
equation [1]. Its derivation starts from a field theoretical identity for the 4-point func- 
tion, the Dyson-Schwinger equation represented in Fig. 2. The kernel is given by the 
sum of all A5-irreducible diagrams, i.e., diagrams that cannot be separated into two 
disconnected parts by cutting one internal A- and one 5-propagator in any possible 
way. The identity becomes obvious after solving it iteratively, which yields an expan- 
sion in powers of the kernel, also represented in Fig. 2. A pole at P 2 = M 2 of the 4-point 
function can then be shown to exist if and only if there exists a solution to the (homo- 
geneous) Bethe-Salpeter equation depicted in Fig. 3 (compare with Figs. 1 and 2), with 



In practice, calculating the kernel exactly is as impossible as calculating the 4- 
point function itself, and the practical value of the Bethe-Salpeter equation consists in 
suggesting an approximation scheme appropriate for bound state calculations, in other 
words, it suggests a choice of a subset of Feynman diagrams to be summed up through an 
approximation of the kernel, and at the same time it devises a method for its summation, 
namely the integration of the equation with the approximate kernel. 

Due to the complexity in solving 4-dimensional integral equations, virtually all equa- 
tions that actually have been solved use the so-called ladder approximation, where the 
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FIGURE 2. An exact relation for the 4-point function, together with its iterative solution in powers of 
the AB-irreducible kernel (represented as a hatched box). 




FIGURE 3. Diagrammatic representation of the Bethe-Salpeter equation. The propagators with mo- 
menta pa and pg have to be understood as amputated, whereas they are contained in y(pa,Pb)- 



kernel is replaced by the one-particle exchange diagram. The iterative solution of the 
Dyson-Sch winger equation for the 4-point function (Fig. 2) in this approximation then 
leads to a sum of ladder diagrams, with the propagators of the exchanged particles as 
rungs. In addition, all self-energy corrections are neglected. In the simplest case of three 
different scalar particles with cubic couplings, the Bethe-Salpeter equation becomes ex- 
plicitly 
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where /i is the mass of the exchanged particle and g the coupling constant. 

Unfortunately, the ladder approximation has a number of well-known defects. To 
begin with, the covariant 4-dimensional formulation implies a dynamics in the relative 
energy p° = p A — p B which is related to the existence of so-called abnormal solutions [3, 
4]. These solutions have no non-relativistic counterparts and appear to be inappropriate 
as solutions of the bound state problem. In the one-body limit nig — > °°, the Bethe- 
Salpeter equation in the ladder approximation does not reduce to a relativistic equation 
for particle A in the field generated by the fixed source B, as is dictated by physical 



considerations [4, 5]. Furthermore, gauge invariance and s-u crossing symmetry are 
violated. 

It is also known that, to remedy the latter three defects, besides the one-particle 
exchange all crossed ladder diagrams have to be included in the kernel, defined as the set 
of all diagrams in which several bosons are exchanged between the two constituents in 
an A5-irreducible way (the simplest example being the crossed box diagram). Needless 
to say, this infinite subset of diagrams cannot be summed exactly. 

Together with the technical difficulties in resolving the 4-dimensional equation, these 
problems have lead to the development of so-called quasi-potential equations where 
the relative energy po is integrated over in such a way that covariance is conserved 
(although not manifestly). While this leads to 3-dimensional equations and eliminates 
the possibility of abnormal solutions, a consideration of the pole structure of the ladder 
and crossed ladder diagrams at the same time allows to approximately sum up all these 
contributions to the kernel. The crucial observation is that the crossed ladder diagrams 
approximately cancel some of the poles of the ladder diagrams, so that an improved 
kernel can be obtained by taking into account less residues in the po-integration [5, 6]. 

Out of the infinity of possible quasi-potential equations, I will cite the three most 
popular choices, all of which are consistent in the one-body limit where some of the pole 
cancellations become exact. The first quasi-potential equation to be formulated was the 
Blancenbecler-Sugar-Logunov-Tavkhelidze (BSLT) equation [7]. It is the "minimal" 
choice in taking two-particle unitarity into account via a dispersion relation. In the c.m.s. 
and for equal constituent masses, the equation reads in the scalar model theory, 

[4(p2+mV£V(p) _, 2 /^^_^L_.0 (4, 

with the lowest-order approximation to the kernel. I have put the equation in a form that 
resembles the non-relativistic Schrodinger equation, and E = M = \fP 2 . 

The next quasi-potential equation to be presented is the Gross (or spectator) equation 
[6]. It is obtained by putting one of the constituents on its mass shell by hand, and reads 
for the scalar model considered before, 
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(p-q) 2 - (vV + ™ 2 - \/q 2 + m 2 ) 2 + /i 2 

The difference of kinetic energies appearing in the denominator of the potential term has 
its origin in the retardation of the interaction through scalar boson exchange. 

Finally, I cite the equal-time equation [8], which results from neglecting the relative 
energy in the kernel of the Bethe-Salpeter equation, and in addition takes the contribu- 
tion of a pole from the crossed box diagram into account. Explicitly, 
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FIGURE 4. The ground state mass mo as a function of the dimensionless coupling constant for jj. jm = 
0.15, and the ground state (equal-time) wave function for iuq = 1.882m; results from the Bethe-Salpeter 
equation (BSE), the BSLT, equal-time (ET) and Gross (with and without retardation) quasi-potential 
equations, and from the Feynman-Sch winger representation (FSR) [10]. 
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As the last method within the Lagrangian approach to bound states I will discuss a 
different technique, inspired by the worldline formulation of QFT rather than by per- 
turbation theory. In this method, the field-theoretical Euclidean functional integral for 
the 4-point function is rewritten as a path integral over the trajectories of the A- and B- 
particles in the so-called Feynman-Sch winger representation. The field corresponding 
to the exchanged boson can then be analytically integrated out [9]. 

In the quenched approximation without closed A- and 5-loops, and after neglecting 
self-energy corrections to the A- and 5-propagators and vertex corrections, the 4-point 
function G(yA,yB',XA,XB) is given by the sum of all ladder and crossed ladder graphs 
(not only A5-irreducible ones). It can be calculated numerically by discretization of the 
path integrals and Monte Carlo methods. The mass Mo of the ground state can then be 
extracted from the dependence of the 4-point function on T = (y^ + — — x$) /2 
in the limit T — > °°, while the dependence on y& — yt? in the same limit determines the 
(relative) Bethe-Salpeter wave function. The numerical results establish, for the first 
time, a benchmark for bound state calculations, at least as long as A- and 5-loops, self 
energies and vertex corrections are neglected. 

The results for the ground state of the scalar model theory with mass ratio ju/m = 
0.15 are depicted in Fig. 4 for the Bethe-Salpeter equation, the BSLT, equal-time 
and Gross equations and from the Feynman-Schwinger representation [10]. Included 
is also a calculation with the Gross equation but neglecting the retardation effects, and 
a representation of the unphysical branch of the solution of the Gross equation (with 
retardation). It is seen that all equations underestimate the binding energy (coming from 
ladder and crossed ladder graphs), where the Bethe-Salpeter equation gives the "worst" 
and the equal-time equation the relatively "best" results. The inclusion of retardation 



and the additional pole from the crossed box improve the results. 

Fig. 4 also compares the different results for the ground state wave functions for a 
fixed ground state mass of mo = 1.882m (and hence varying coupling constants). As far 
as the quality of the approximations is concerned, the outcome is precisely opposite to 
the comparison of the binding energies. In conclusion, there is no compelling evidence 
to prefer one of the quasi-potential equations over the others. 

Turning now to the Hamiltonian approach to QFT, a bound state (or any scattering 
state, for that matter) is defined to be an eigenstate of the Hamiltonian H. The Hamilto- 
nian is decomposed into a free part Hq and an interaction part Hj, where the eigenstates 
of Ho are well-known and used to construct the Fock space & of the theory. It is as- 
sumed that any state in the interacting theory has a unique expansion in terms of the 
Fock space states. Then the fundamental equation to be solved is 

H\ ¥ )=E\ ¥ ), (7) 

Now in general Hj mixes the different n-particle sectors of so that Eq. (7) is really 
an infinite system of coupled (integral) equations for all sectors of Fock space, and 
the eigenstates contain contributions from every sector (except where forbidden by 
conservation laws). A systematic approximation scheme for the solution of Eq. (7) is 
hence called for. 

Most of the calculations in a Hamiltonian formalism are currently done in QFT 
quantized on the light front (the 3-dimensional space perpendicular to an arbitrarily 
chosen direction on the forward light cone). The light-front dynamics differs from the 
usual equal-time quantized theory in several aspects: the "spectrum condition" p + > 
which is related to the form of the light-front Hamiltonian, holds for real and virtual 
particles, and together with p + -conservation at the vertices implies that there are no 
pair production processes in the vacuum. However, the vacuum is, contrary to earlier 
folklore, in general not trivial due to zero modes (constant field configurations on the 
light front). In particular, this seems to be the case in non-abelian gauge theories [11]. 
Still, the structure of the vacuum should be much simpler than in the usual formulation. 
Of equal importance is the suppression of higher Fock components by kinematical 
factors, even in the case of rather large coupling constants. Obviously, these properties 
are particularly relevant to QCD. 

A serious problem is posed by the longitudinal divergencies which require non-local 
counterterms, i.e., entire functions of p + instead of a few numbers. A new, "non- 
perturbative", renormalization theory has to be developed, possibly along the lines of 
the coupling coherence scheme proposed by Perry and Wilson [12]. 

The central idea behind all the approaches to the (approximate) solution of Eq. (7) in 
light-front quantization is the "reduction of Fock space" by similarity transformations, 
i.e., the block diagonalization of the Hamiltonian converting (7) into an (approximately) 
equivalent equation on a small subspace of the Fock space & . In practice, the similarity 
transformations are usually realized as a sequence of unitary transformations. Several 
similarity transformation schemes have been developed [13]. Among the few applica- 
tions of this technique to given field theories I mention the calculation of the fine and 
hyperfine structure of positronium for a larger fine structure constant a = 0.3 in the dif- 
ferent schemes [14]. The numerical results are satisfactory, however, a (quantitatively 



small) violation of rotational invariance and a dependence on the similarity transfor- 
mation scheme remain. The situation is much less clear for bound state calculations in 
Yukawa theory because of renormalization problems [11, 15]. 

Given the difficulties with theories quantized on the light-front, one might wonder 
whether the similarity transformation method for the solution of Eq. (7) could not 
be applied successfully to the usual equal-time quantized theory. In fact, a related 
formalism is currently being developed. It starts from the following generalization of the 
Gell-Mann-Low theorem in QFT [2] : the time evolution determined by the Hamiltonian 
H(t) = Hq + e~ £ ^Hi in the adiabatic limit £ — > 0, maps eigenstates of Hq at t — > — °° 
that are isolated in the spectrum to eigenstates of H at t = 0. The same time evolution 
operator can be shown to map //o-invariant linear subspaces of & to //-invariant 
subspaces under more general conditions [16]. The Hamiltonian can then be similarity 
transformed to the //o-invariant subspace. For example, if one considers the subspace of 
all free 2-particle states, then Eq. (7) is equivalent to a quantum-mechanical 2-particle 
Schrodinger equation (on the subspace). 

Applied to the scalar model considered before, the similarity transformation leads to 
second order in g to the following effective Schrodinger equation for A5-bound states 
[17]: 
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In the hermitian version that was used by Wilson [18] and by Glockle and Muller [19], 
the equation becomes 
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A comparison of the numerical results for the binding energy of the ground state 
from the solution of Eq. (8) with those of the Bethe-Salpeter equation and several other 
equations which have not been discussed here for different reasons, is presented in Fig. 5 
for the case fx = [17]. Solutions of Eq. (8) for some excited s-states are also shown. The 
ladder approximation of the Bethe-Salpeter equation is known in this case as the Wick- 
Cutkosky model, and is exactly solvable [4]. The results from Eq. (8) lie roughly half- 
way between the Bethe-Salpeter and the non-relativistic results (value (—1) in Fig. 5), 
which is reasonable in view of the results from the Fey nman-Sch winger representation 
shown in Fig. 4. 
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FIGURE 5. The binding energy = M — 2m of the ground state and several excited s-states as a 
function of the "fine structure constant" a = g 2 /{\6nm 2 ) for jU = 0, normalized to the non-relativistic 
ground state energy [17] (M r = m/2 is the reduced mass); results from the effective Schrodinger equation 
(8) ("Bloch- Wilson"), the ladder approximation of the Bethe-Salpeter equation ("Wick-Cutkosky"), a 
straightforward Fock space expansion [20], the leading Regge trajectory [21], the Bethe-Salpeter equation 
in covariant light-front dynamics [22], and the Haag expansion [23]. 
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